Moderate Resolution Imaging Spectroradiometer Research Articles

Accelerated reduction in China's cropland fires against the background of policy enhancement

Cropland fires contribute significantly to the emission of air pollutants, with implications for the atmospheric environment, public health, and climate change. However, the development of cropland fires in response to straw burning policy interventions has not been fully explored. With the Moderate Resolution Imaging Spectroradiometer Collection 6 (MODIS C6) active fire products and China land-cover data (CLCD), counts and fire radiative power (FRP) of cropland fire spots were extracted in China from 2003 to 2022, followed by analyses of their spatiotemporal patterns and evolutionary characteristics at the provincial, regional, national, and grid (0.25° × 0.25°) levels and at monthly to annual scales via trend analysis, kernel density analysis, and geographic distribution analysis. The key findings include: (1) National policies banning straw burning reversed the changes in the quantity and intensity of China's cropland fires in 2013, with the frequency and FRP in 2022 declining by 58% and 62% compared with the 2014 peak. (2) Cropland fires were highly prevalent in Northeast and North China (NEC and NC), with significant regional differences, particularly in spring, summer and autumn. (3) A notable lag effect of burning ban policies was seen in Northeast China (NEC), where the increasing trends of the frequency and intensity of cropland fires were not reversed until 2017, followed by declines in frequency and FRP of 67% and 62%, respectively during 2017–2022. (4) Cropland fires showed an increasing concentration trend but with enormous regional differences, with an increasing trend of a “high-high” concentration pattern in NEC and a decreasing trend in NC. This study contributes to understanding the effectiveness of national policies that ban straw burning and support cropland fire management and policy-making.

Development of long‐term spatiotemporal continuous NDVI products for alpine grassland from 1982 to 2020 in the Qinghai–Tibet Plateau, China

AbstractBackgroundThe time‐series data of the Normalized Difference Vegetation Index (NDVI) is a crucial indicator for global and regional vegetation monitoring. However, the current assessment of global and regional long‐term vegetation changes is subject to large uncertainties due to the lack of spatiotemporally continuous time‐series data sets.MethodsIn this study, a long time‐series monthly NDVI data set with a spatial resolution of 250 m from 1982 to 2020 was developed by combining Moderate Resolution Imaging Spectroradiometer (MODIS) and AVHRR (Advanced Very High‐Resolution Radiometer) time‐series NDVI products using the Random Forest (RF) downscaling model.ResultsCompared to the MODIS NDVI product, the fused product shows RMSE and mean absolute error ranging from 0 to 0.075 and from 0 to 0.05, respectively, with R2 values mostly above 0.7.ConclusionsThe long time‐series NDVI products generated in this study are reliable in terms of accuracy and have great potential for long‐term dynamic monitoring of terrestrial ecosystems on the Qinghai–Tibet Plateau.

Relationship between aerosol and cloud characteristics over Delhi in North India during the dry and wet season

The study investigates the spatio-temporal variations of aerosol and cloud properties retrieved from Moderate Resolution Imaging Spectroradiometer (MODIS) satellite at an urban megacity Delhi in the northern part of India during the dry and wet season for 15 years period from November 2004 to October 2019. This research comprises the study of crucial aerosol optical properties such as aerosol optical depth (AOD), Angstrom exponent (AE) along with the cloud properties such as water vapor (WV), cloud liquid water path (CLWP), cloud effective radius liquid (CERL), cloud fraction (CF), cloud optical depth (COD), cloud top pressure (CTP), and cloud top temperature (CTT). High AOD was observed across the Indo-Gangetic Plain (IGP), including Delhi and exhibited a range of 0.49–0.89 during the dry season and 0.73 to 0.98 during the wet season with a mean value of 0.73 ± 0.05 during the entire study period. High AE value (mean 1.40 ± 0.05) indicates dominance of fine-mode aerosols over the station, associated with high anthropogenic activities. A significantly large increasing trend in AOD was observed over the station during the wet season as compared to the dry season. The seasonal studies represent the maximum value of WV and CTP in the wet season. In addition, CF, CERL, COD and CLWP values were also high during the wet season. Furthermore, a correlation analyses between the aerosol and cloud parameters were carried out during both the seasons. From the analysis, AOD shows a positive correlation with WV, CF, COD and CLWP in both the seasons whereas a negative correlation was observed with CTP in wet season, while positive in dry season. On the other hand, a poor correlation was observed between AOD verses AE and CTT in both the seasons. These findings can contribute to advancing our understanding of aerosol-cloud interactions and inflate air quality assessments at an urban megacity like Delhi.

Temporal and spatial variations in dust activity in Australia based on remote sensing and reanalysis datasets

Abstract. Spatial and temporal variations in the level of dust activity can provide valuable information for policymaking and climate research. Recently, Moderate Resolution Imaging Spectroradiometer (MODIS) aerosol products have been successfully used for retrieving dust aerosol optical depth (DAOD), especially over bright dust source areas, and Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) aerosol reanalysis provides DAOD and additionally other dust-aerosol-related parameters. In this study, spatial and temporal variations in dust activity in Australia were analyzed using MODIS and MERRA-2 combined (M&M) DAOD and MERRA-2 near-surface dust concentrations and estimated PM10 for the period 1980–2020. Validation results show that M&M DAOD has an expected error of ±(0.016+0.15τ) compared to the ground observations at the AErosol RObotic NETwork (AERONET) dust storms on populated areas sites. MERRA-2 near-surface dust concentrations show a power-law relationship with visibility data collected at meteorological stations with an r2 value from 0.18 to 0.44, and the estimated MERRA-2 PM10 shows similar temporal variations and correlates with ground-based PM10 data with an r2 value from 0.14 to 0.44 at six selected stations in Australia. Moreover, MERRA-2 horizontal dust flux shows the same major dust pathways as those in previous studies and similar dust emissions and deposition areas identified using ground-based observations. Dust events based on DAOD over eastern Australia are concentrated in the north in December, are concentrated in the south in February, and can occur anywhere in January. Near-surface dust concentration was found to be the highest (over 200 µg m−3) over the center of Lake Eyre basin in central Australia and radially decreased toward the coast to below 20 µg m−3 via the two main pathways in the southwest and northeast. The ratio of near-surface dust concentration to PM10 shows a similar spatial pattern. Total dust emission was estimated to be 40 Mt (megatonnes) per year over the period 1980–2020, of which nearly 50 % was deposited on land and the rest exported away from the Australian continent.

Introduction to the NJIAS Himawari-8/9 Cloud Feature Dataset for climate and typhoon research

Abstract. The use of remote sensing methods to accurately measure cloud properties and their spatiotemporal changes has been widely welcomed in many fields of atmospheric research. The Nanjing Joint Institute for Atmospheric Sciences (NJIAS) Himawari-8/9 Cloud Feature Dataset (HCFD) provides a comprehensive description of cloud features over the East Asia and west North Pacific regions for the 7-year period from April 2016 to December 2022. Multiple cloud variables, such as cloud mask, phase/type, top height, optical thickness, and particle effective radius, as well as snow, dust, and haze masks, were generated from the visible and infrared measurements of the Advanced Himawari Imager (AHI) on board the Japanese geostationary satellites Himawari-8 and Himawari-9 using a series of recently developed cloud retrieval algorithms. Verifications with the Cloud–Aerosol Lidar with Orthogonal Polarization (CALIOP) 1 km cloud layer product and the Moderate Resolution Imaging Spectroradiometer (MODIS) Level-2 cloud product (MYD06) demonstrate that the NJIAS HCFD gives higher skill scores than the Japanese Himawari-8/9 operational cloud product for all cloud variables except for cloud particle effective radius. The NJIAS HCFD even outperforms the MYD06 in nighttime cloud detection; cloud-top height, pressure, and temperature estimation; and infrared-only cloud-top phase determination. All evaluations are performed at the nominal 2 km scale, not including the effects of sub-pixel cloudiness or very thin cirrus. Two examples are presented to demonstrate applications of the NJIAS HCFD for climate and typhoon research. The NJIAS HCFD has been published in the Science Data Bank (https://doi.org/10.57760/sciencedb.09950, Zhuge 2023a; https://doi.org/10.57760/sciencedb.09953, Zhuge 2023b; https://doi.org/10.57760/sciencedb.09954, Zhuge 2023c; https://doi.org/10.57760/sciencedb.10158, Zhuge 2023d; https://doi.org/10.57760/sciencedb.09945, Zhuge 2023e).

The Effect of Physically Based Ice Radiative Processes on Greenland Ice Sheet Albedo and Surface Mass Balance in E3SM

AbstractA significant portion of surface melt on the Greenland Ice Sheet (GrIS) is due to dark ice regions in the ablation zone, where solar absorption is influenced by the physical properties of the ice, light absorbing constituents (LACs), and the overlying crustal surface or melt ponds. Earth system models (ESMs) typically prescribe the albedo of ice surfaces as a constant value in the visible and near‐infrared spectral regions. This work advances ESM ice radiative transfer modeling by (a) incorporating a physically based radiative transfer model (SNow, ICe and Aerosol Radiation model Adding‐Doubling Version 4; SNICAR‐ADv4) into the Energy Exascale Earth System Model (E3SM), (b) determining spatially and temporally varying bare ice physical properties over the GrIS ablation zone from satellite observations to inform SNICAR‐ADv4, and (c) assessing the impacts on simulated GrIS albedo and surface mass balance associated with modeling of more realistic bare ice albedo. GrIS‐wide bare ice albedo in E3SMv2 is overestimated by ∼4% in the visible and ∼7% in the near‐infrared wavelengths compared to the Moderate Resolution Imaging Spectroradiometer. Our bare ice physical property retrieval method found that LACs, ice crustal surfaces, and melt ponds reduce visible albedo by 30% in the bare ice region of the GrIS ablation zone. The realistic bare ice albedo reduces surface mass balance by ∼145 Gt, or 0.4 mm of sea‐level equivalent between 2000 and 2021 compared to the default E3SM. This work highlights the importance of simulating bare ice albedo accurately and realistically to improve our ability to quantify changes in the GrIS surface mass and radiative energy budgets.

Interannual Variations in the Summer Coastal Upwelling in the Northeastern South China Sea

This study scrutinizes interannual (2003–2023) variations in coastal upwelling along the Guangdong Province during summers (June–August) in the northeastern South China Sea (NESCS) by comprehensively applying the moderate-resolution imaging spectroradiometer (MODIS) remote sensing sea surface temperature (SST) and chlorophyll concentration (CHL) data and the model reanalysis product. The results show that SST and upwelling intensity in the sea area have significant (p < 0.05) rising trends in the last 21 years. The CHL shows an upward but insignificant trend, which is affected simultaneously by the rise in SST and the enhancement of upwelling. Further analysis reveals that the interannual variations in upwelling are robustly related to the wind fields’ variations in the coastal region. A clockwise/counter-clockwise anomaly in the wind field centered on the NESCS facilitates alongshore/onshore winds near the Guangdong coast, which can strengthen/weaken coastal upwelling. Based on the correlation between wind field variations and large-scale climate factors, long-term variations in the upwelling intensity can be primarily predicted by the Oceanic Niño Index.

Parameterizing spectral surface reflectance relationships for the Dark Target aerosol algorithm applied to a geostationary imager

Abstract. Originally developed for the moderate resolution imaging spectroradiometer (MODIS) in polar, sun-synchronous low earth orbit (LEO), the Dark Target (DT) aerosol retrieval algorithm relies on the assumption of a surface reflectance parameterization (SRP) over land surfaces. Specifically for vegetated and dark-soiled surfaces, values of surface reflectance in blue and red visible-wavelength bands are assumed to be nearly linearly related to each other and to the value in a shortwave infrared (SWIR) wavelength band. This SRP also includes dependencies on scattering angle and a normalized difference vegetation index computed from two SWIR bands (NDVISWIR). As the DT retrieval algorithm is being ported to new sensors to continue and expand the aerosol data record, we assess whether the MODIS-assumed SRP can be used for these sensors. Here, we specifically assess SRP for the Advanced Baseline Imager (ABI) aboard the Geostationary Operational Environmental Satellite (GOES)-16/East (ABIE). First, we find that using MODIS-based SRP leads to higher biases and artificial diurnal signatures in aerosol optical depth (AOD) retrievals from ABIE. The primary reason appears to be that the geostationary orbit (GEO) encounters an entirely different set of observation geometry than does LEO, primarily with regard to solar angles coupled with fixed-view angles. Therefore, we have developed a new SRP for GEO that draws the angular shape of the surface bidirectional reflectance. We also introduce modifications to the parameterization of both red–SWIR and blue–red spectral relationships to include additional information. The revised red–SWIR SRP includes the solar zenith angle, NDVISWIR, and land-type percentage from an ancillary database. The blue–red SRP adds dependencies on the scattering angle and NDVISWIR. The new SRPs improve the AOD retrieval of ABIE in terms of overall less bias and mitigation of the overestimation around local noon. The average bias of the DT AOD compared to the Aerosol Robotic Network (AERONET) AOD shows a reduction from 0.08 to 0.03, while the bias of local solar noon decreases from 0.12 to 0.03. The agreement between the DT and AERONET AOD is established through a regression slope of 1.06 and a y intercept of 0.01 with a correlation coefficient of 0.74. By using the new SRP, the percentage of data falling within the expected error range (±0.05 % + 15 %) is notably increased from 54 % to 78 %.

Seasonal variations in chlorophyll–a and sea surface temperature in the exclusive economic zone of Sri Lanka

Oceanic phytoplankton plays a crucial role in regulating the biogeochemical cycle in marine ecosystems. Light harvesting pigments such as chlorophyll–a (chl–a) can be used to understand marine primary productivity, especially phytoplankton abundance. Sea surface temperature is closely linked with phytoplankton abundance and climatic changes. The objective of the current study is to determine spatial and temporal variations of chl–a and sea surface temperature and their co-variation, over 12 years using statistical methodologies and Remote Sensing, Geographic Information System (GIS) and Cloud Computing Approach. In this study, the available sea surface temperature and chl–a data acquired from a Moderate Resolution Imaging Spectroradiometer (MODIS) satellite sensor were proceeded through Google Earth Engine and Earth Science Data System. Data visualization was carried out using ArcMap™ software. Results show a weak negative linear relationship (–0.40, P < 0.05) between sea surface temperature and chl–a. Therefore, additional favorable environmental factors such as coastal upwelling, and nutrient discharge through river run-off can control the primary productivity of phytoplankton. We also illustrated spatial and temporal maps for chl–a and sea surface temperature distribution in the exclusive economic zone of Sri Lanka. Chlorophyll–a concentration is enhanced during the southwest monsoon, and the lowest concentration is observed during the first inter-monsoon period. Next, the highest sea surface temperature is recorded during the first inter-monsoon period, and the lowest sea surface temperature is observed during the northeast monsoon. Seasonal behavior of chl–a concentration and sea surface temperature are changed as chl–a (SWM) > chl–a (NEM) > chl–a (SIM) > chl–a (FIM) and SST (FIM) > SST (SIM) > SST (SWM) > SST (NEM), respectively. In addition, both trend lines for each variable show a positive gradient, which may be attributed to global warming. Finally, we observed the alteration of chl–a and sea surface temperature that linked to changes in ocean chemistry during the global COVID–19 pandemic.

Modelling Floodplain Vegetation Response to Climate Change, Using the Soil and Water Assessment Tool (SWAT) Model Simulated LAI, Applying Different GCM’s Future Climate Data and MODIS LAI Data

Scientists widely agree that anthropogenically driven climate change significantly impacts vegetation growth, particularly in floodplain areas, by altering river flow and flood regimes. This impact will accelerate in the future, according to climate change projections. For example, in Australia, climate change has been attributed to a decrease in winter precipitation in the range of 56% to 72.9% and an increase in summer from 11% to 27%, according to different climate scenarios. This research attempts to understand vegetation responses to climate change variability at the floodplain level. Further, this study is an effort to enlighten our understanding of temporal climate change impacts under different climate scenarios. To achieve these aims, a semi-distributed hydrological model was applied at a sub-catchment level to simulate the Leaf Area Index (LAI). The model was simulated against future time series of climate data according to Global Climate Model (GCM) projections. The time series data underwent a non-parametric Mann–Kendall test to detect trends and assess the magnitude of change. To quantify the model’s performance, calibration and validation were conducted against the Moderate Resolution Imaging Spectroradiometer (MODIS) LAI. The calibration and validation results show Nash–Sutcliffe efficiency (NSE) values of 0.85 and 0.78, respectively, suggesting the model’s performance is very good. The modeling results reveal that the rainfall pattern fluctuates under climate projections within the study site, in which vegetation tends to be more vibrant during the warmer seasons. Moreover, the modeling results highlighted increases in the average projected future winter temperatures, which can help vegetation growth during winter. The results of this study may be employed for sustainable floodplain management, restoration, land-use planning, and policymaking, and help floodplain communities better prepare for and respond to changing flood patterns and related challenges under a future changing climate.

Retrieval of sea ice drift in the Fram Strait based on data from Chinese satellite HaiYang (HY-1D)

Abstract. Melting of sea ice in the Arctic has accelerated due to global warming. The Fram Strait (FS) serves as a crucial pathway for sea ice export from the Arctic to the North Atlantic Ocean. Monitoring sea ice drift (SID) in the FS provides insight into how Arctic sea ice responds to the climate change. The SID has been retrieved from Sentinel-1 synthetic aperture radar (SAR), Advanced Very High Resolution Radiometer (AVHRR), Moderate Resolution Imaging Spectroradiometer (MODIS), and Advanced Microwave Scanning Radiometer for EOS (AMSR-E), and further exploration is needed for the retrieval of SID using optical imagery. In this paper, we retrieve SID in the FS using the Chinese HaiYang1-D (HY-1D) satellite equipped with the Coastal Zone Imager (CZI). A multi-template matching technique is employed to calculate the cross-correlation, and subpixel estimation is used to locate displacement vectors from the cross-correlation matrix. The dataset covering March to May 2021 was divided into hourly and daily intervals for analysis, and validation was performed using Copernicus Marine Environment Monitoring Service (CMEMS) SAR-based product and International Arctic Buoy Programme (IABP) buoy. A comparison with the CMEMS SID product revealed a high correlation with the daily interval dataset; however, due to the spatial and temporal variability in the sea ice motion, differences are observed with the hourly interval dataset. Additionally, validation with the IABP buoys yielded a velocity bias of −0.005 m s−1 and RMSE of 0.031 m s−1 for the daily interval dataset, along with a flow direction bias of 0.002 rad and RMSE of 0.009 rad, respectively. For the hourly interval dataset, the velocity bias was negligible (0 m s−1), with a RMSE of 0.036 m s−1, while the flow direction bias was 0.003 rad, with a RMSE of 0.010 rad. In addition, during the validation with buoys, we found that the accuracy of retrieving the SID flow direction is distinctly interrelated with the sea ice displacement.

Trends and drivers of Arctic-boreal fire intensity between 2003 and 2022

Climate change has disproportional effects on Arctic-boreal ecosystems, as the increase of air temperatures in these northern regions is several times higher than the global average. Ongoing warming and drying have resulted in recent record-breaking fire years in Arctic-boreal ecosystems, resulting in substantial carbon emissions that might accelerate climate change. While recent trends in Arctic-boreal burned area have been well documented, it is still unclear how fire intensity has changed. Fire intensity relates to the energy release from combustion and to a large extent drives the impact of a fire on the vegetation and soils, the emission of various gasses and the combustion completeness of different fuels. Here, we used the active fire product from the Moderate Resolution Imaging Spectroradiometer (MODIS) to examine trends in fire radiative power (FRP) over the entire Arctic-boreal region. We found a significant increase in annual median fire intensity between 2003 and 2022 in Eurasian boreal forests, for both daytime (increase of 0.392 MW/km2 per year, R2 = 0.56, p < 0.001) and nighttime fires (increase of 0.175 MW/km2 per year, R2 = 0.47, p < 0.001), while no general trend in FRP was observed in boreal North America. This increase in FRP in Eurasian boreal forests was strongly associated with simultaneous increases in air temperature, vapour pressure deficit, fire weather and fuel availability. We estimated that for Eurasia with each degree increase in air temperature, annual median daytime FRP increases with 1.58 MW/km2 in the tundra and 0.94 MW/km2 in the taiga. Climate change has thus resulted in a widespread and clear increase in fire intensity in central and eastern Eurasia while we could not discern clear trends in Arctic-boreal North America. Arctic-boreal fire intensity may further increase with climate change, with potentially major consequences for fire regimes, carbon emissions and society.

Arctic Greening Trends: Change Points in Satellite-Derived Normalized Difference Vegetation Indexes and Their Correlation with Climate Variables over the Last Two Decades

In this study, we utilized NDVI data from the moderate resolution imaging spectroradiometer (MODIS) alongside climatic variables obtained from a reanalyzed dataset to analyze Arctic greening during the summer months (June–September) of the last two decades. This investigation entailed a detailed analysis of these changes across various temporal scales. The data indicated a continuous trend of Arctic greening, evidenced by a 1.8% per decade increment in the NDVI. Notably, significant change points were identified in June 2012 and September 2013. A comparative assessment of NDVI pre- and post-these inflection points revealed an elongation of the Arctic greening trend. Furthermore, an anomalous increase in NDVI of 2% per decade was observed, suggesting an acceleration in greening. A comprehensive analysis was conducted to decipher the correlation between NDVI, temperature, and energy budget parameters to elucidate the underlying causes of these change points. Although the correlation between these variables was relatively low throughout the summer months, a distinct pattern emerged when these periods were dissected and examined in the context of the identified change points. Preceding the change point, a strong correlation (approximately 0.6) was observed between all variables; however, this correlation significantly diminished after the change point, dropping to less than half. This shift implies an introduction of additional external factors influencing the Arctic greening trend after the change point. Our findings provide foundational data for estimating the tipping point in Arctic terrestrial ecosystems. This is achieved by integrating the observed NDVI change points with their relationship with climatic variables, which are essential in comprehensively understanding the dynamics of Arctic climate change, particularly with alterations in tundra vegetation.

HiQ-LAI: a high-quality reprocessed MODIS leaf area index dataset with better spatiotemporal consistency from 2000 to 2022

Abstract. Leaf area index (LAI) is a crucial parameter for characterizing vegetation canopy structure and energy absorption capacity. The Moderate Resolution Imaging Spectroradiometer (MODIS) LAI has played a significant role in landmark studies due to its clear theoretical basis, extensive historical time series, extensive validation results, and open accessibility. However, MODIS LAI retrievals are calculated independently for each pixel and a specific day, resulting in high noise levels in the time series and limiting its applications in the regions of optical remote sensing. Reprocessing MODIS LAI predominantly relies on temporal information to achieve smoother LAI profiles with little use of spatial information and may easily ignore genuine LAI anomalies. To address these problems, we designed the spatiotemporal information compositing algorithm (STICA) for the reprocessing of MODIS LAI products. This method integrates information from multiple dimensions, including pixel quality information, spatiotemporal correlation, and the original retrieval, thereby enabling both “reprocessing” and “value-added data” with respect to the existing MODIS LAI products, leading to the development of the high-quality LAI (HiQ-LAI) dataset. Compared with ground measurements, HiQ-LAI shows better performance than the original MODIS product with a root-mean-square error (RMSE) or bias decrease from 0.87 or −0.17 to 0.78 or −0.06, respectively. This is due to the improvement of HiQ-LAI with respect to capturing the seasonality in vegetation phenology and reducing abnormal time-series fluctuations. The time-series stability (TSS) index, which represents temporal stability, indicated that the area with smooth LAI time series expanded from 31.8 % (MODIS) to 78.8 % (HiQ) globally, and this improvement is more obvious in equatorial regions where optical remote sensing cannot usually achieve good performance. We found that HiQ-LAI demonstrates superior continuity and consistency compared with raw MODIS LAI from both spatial and temporal perspectives. We anticipate that the global HiQ-LAI time series, generated using the STICA procedure on the Google Earth Engine (GEE) platform, will substantially enhance support for diverse global LAI time-series applications. The 5 km 8 d HiQ-LAI datasets from 2000 to 2022 are available at https://doi.org/10.5281/zenodo.8296768 (Yan et al., 2023).

Three-dimensional reconstruction of aerosols based on spectral radiance matching

This paper evaluates a three-dimensional scene construction algorithm that assigns CALIPSO-observed profiles to receiving pixels in MODIS imagery through spectral radiance matching (SRM), so that a limited number of vertically resolved aerosol observations can be expanded into a regionally continuous three-dimensional aerosol structure. The SRM algorithm was tested using moderate resolution imaging spectro-radiometer (MODIS) and cloud-aerosol lidar and infrared pathfinder satellite observation (CALIPSO) data for the period of April 10–25, 2015. First, the spectral radiance channel and the surface reflectance channel were selected for combination tests to select the best channel combination 01-04-07-26 for the application of SRM algorithm in the field of aerosol three-dimensional reconstruction. The average aerosol matching rate for this channel combination can reach 90.57%, 81.64%, 78.6%, 74.12%, and 65.96% at extended distances of 5, 20, 30, 50, and 100 km, respectively. In addition, the combination of the average correct match rate and the reconstructed aerosol structure revealed that the most suitable extension distance is up to 100 km. Then, for the reconstruction profile validation, lidar sites under four weather conditions (clean sky, aerosol loading, thin clouds, and thick clouds) were selected for validation. The comparison results reveal that the average relative errors are about 30%, 60%, 30%–60%, and 70% for these four conditions, respectively. Overall, the average relative error of the reconstructed profiles within 100 km of the orbital track of CALIPSO (except under thick cloud conditions) is less than 60%, which is within a reasonable error envelop. There are also some limitations that can be improved later by adding more constraints: 1) it is difficult to reconstruct some short-time, small-scale aerosol or cloud features; 2) it is difficult to reconstruct some hierarchical structures under thick clouds due to signal attenuation caused by thick clouds; 3) there is an underestimation of the attenuated backscatter coefficient under clouds. Finally, for aerosol optical depth (AOD), the SRM algorithm largely underestimates AOD. The reconstructed AOD is in good agreement with aerosol robotic network (AERONET), with R2=0.8231 after removal of outliers.

A novel algorithm for ocean chlorophyll-a concentration using MODIS Aqua data

The ability to infer ocean chlorophyll-a concentrations (Chla) from spaceborne instruments is key to assessments of global ocean productivity and monitoring of water quality. Here, we present a novel parametric algorithm, OCG, trained on a set of global in situ high-performance liquid chromatography (HPLC) data that leverages Level-3 remote sensing reflectance (Rrs) products from the Moderate Resolution Imaging Spectroradiometer (MODIS) aboard the Aqua satellite. The OCG algorithm leverages more bands than existing algorithms and also provides pixel-wise uncertainty assessments that enable the calculation of the probability of exceeding specific Chla thresholds. This feature has significant implications for water quality management, particularly in monitoring harmful algal blooms. The OCG surpasses existing algorithms in bias and accuracy without overfitting, especially in coastal areas, where it outperforms the current standard product (CI OC3) by 20 % in median symmetric accuracy. Moreover, the OCG reduces the signed symmetric percentage bias (SSPB) in coastal regions from 41 % (CI OC3) to below 5 %. Globally, the OCG algorithm yields lower Chla in coastal regions, the Southern Ocean and the Mediterranean Sea, and higher values in the open ocean, particularly in ocean gyres and polar regions. For the Chesapeake Bay and the Baltic Sea, for example, daily OCG estimates for 2002 to 2021 are, on average, 2.9 g/L and 3.7 g/L lower than CI OC3 estimates, respectively. The presented approach also shows great potential for other existing and upcoming sensors, enabling widespread application in remote sensing.

Investigating the Role of the Low-Level Jet in Two Winters Severe Dust Rising in Southwest Iran

The dust storms with local and non-local dust sources mostly affect Khuzestan province in southwest (SW) Iran. In this study, the role of the low-level jet in the activation of the internal dust events in SW Iran during two severe dust cases was investigated. For this purpose, the fifth-generation ECMWF reanalysis for the global climate and weather (ERA5) data was used to identify the synoptic patterns and the low-level jet (LLJ) characteristics in the study area. Furthermore, the images of the moderate resolution imaging spectroradiometer (MODIS) sensor, the outputs of the hybrid single-particle Lagrangian integrated trajectory (HYSPLIT) model, and a weather research and forecasting model coupled with chemistry (WRF-Chem) were used to investigate the propagation and transport of the dust particles. The results of the synoptic analysis in both dust cases show the simultaneous occurrence of the divergence zone associated with cyclonic curvature in the subtropical jet stream (STJ) at 300 hPa, causing convergence at 925 hPa, upward motion, and the development of low surface pressure in SW Iran. Examining the vertical wind profile shows the existence of the maximum horizontal wind speeds of 975 to 875 hPa, along with the positive and negative shear below and above it, respectively, which emphasizes the existence of the LLJ and its role in local dust emission. The results of the comparison between the satellite images, WRF-Chem, and HYSPLIT model outputs show the formation and transportation of dust particles from the inner regions of Khuzestan in SW Iran. The horizontal dust surface distribution, vertical raised dust mass, and kinetic energy transfers are well simulated by the WRF-Chem model when LLJ broke at 09:00 to 12:00 UTC. The most important finding of this research is that, for the first time, the role of low-level jet is investigated in the activation of internal dust events in SW Iran.

Assessing topographic effects on forest responses to drought with multiple seasonal metrics from Sentinel-2

Topography determines run-off direction, redistributes groundwater, and affects land surface solar radiation loads and the associated evaporative forcing, consequently, topography can modulate the impact of drought and heat waves on ecosystems. This topographic modulation effect, which typically occurs at the local scale, is often overlooked when assessing ecosystem drought responses using moderate-to-coarse spatial resolution satellite observations, such as the Moderate-resolution Imaging Spectroradiometer (MODIS) imagery. Sentinel-2 and Landsat imagery with finer resolution are suitable for monitoring changes at the local scale, however, studies relying on single vegetation metrics may fail to get a holistic picture of vegetation drought responses, particularly for forests that have complex physiological mechanisms. Here, we performed a comprehensive assessment of the topographic effects on coniferous forest responses to the severe 2018 drought in Scandinavia, using 6 vegetation seasonal metrics during 2017–2021 from the Sentienl-2 High-Resolution Vegetation Phenology and Productivity (HR-VPP) products. We found significant differences (p < 0.05) between sunny and shady aspects, between higher and lower elevations, and between steep and gentle slopes, regarding the maximum impact time, forest drought resistance, and resilience. Specifically, the sunny aspects and steep slopes were related to higher risks of delayed impacts and low resistance and resilience, and elevation and slope were more powerful in regulating the phenology shift and greening rate loss. We also identified different sensitivity in greenness and productivity to topographic effect and greater sensitivity of spring phenology to topographic differences as compared to autumn phenology. The study demonstrates vegetation drought responses represented by multiple seasonal metrics, reveals the prevalent topographic effects at the local scale, and quantifies the magnitudes of the effects with regional statistics.

Spatial and temporal patterns of land surface temperature in Greenland from 2000-2019

Temperature dynamics on the island of Greenland are an important factor in shaping ecological events. Investigating the land surface temperature (LST) patterns is critical for understanding ecological dynamics across different regions. Further melting of the Greenland ice sheet could deva state marine and terrestrial ecosystems. This study used data from Moderate Resolution Imaging Spectroradiometer satellites to understand the seasonal patterns and patterns of LST over the entire island. Focusing on the period between 2000 and 2019, this study used a natural cubic spline model to identify seasonal patterns for all sub-regions. The data were seasonally adjusted and filtered with a second-order autocorrelation component. The spline was fitted again to identify the LST pattern, and a multivariate regression model was then used to adjust for spatial correlation. We illustrate that most of the land surface of Greenland hasstable temperature trends. These observed patterns in LST in Greenland during the study period suggest that the observed ice-sheet melting in Greenland within the last two decades could be due to other factors, not necessarily LST patterns.

Unleashing the power of machine learning and remote sensing for robust seasonal drought monitoring: A stacking ensemble approach

Droughts cause significant economic losses in many regions around the world, highlighting a need to more accurately quantify implications of drought on production and water management. Remote sensing technologies and machine learning-based models offer promising solutions for timely and accurate regional drought monitoring, but the accuracy of such approaches is constrained by both the conceptual design and algorithms underpinning such approaches. We developed a machine learning stacking ensemble approach to overcome such limitations, drawing upon Precipitation Estimation from Remotely Sensed Information using Artificial Neural Network-Climate Data Record (PERSIANN-CDR), Moderate Resolution Imaging Spectroradiometer (MODIS) remote sensing products, and climate zoning data to estimate 3-month scale SPEI (Standardized Precipitation Evapotranspiration Index, SPEI-3) in 9 sub-regions of China. We compared 19 individual machine learning models and used the stacking approach to select the most robust meta-model for SPEI-3 prediction. We found that CatBoost Regressor (CBR), Extra Trees Regressor (ETR), Extreme Gradient Boosting (XGB), Light Gradient Boosting Machine (LGBM), and Random Forest (RF) were the top individual models for predicting drought, while using CBR as the stacked meta-model achieved the best performance. The R2 values for the stacking model with CBR as the meta-model were 0.9065 and 0.8218 in the eastern and western regions, respectively. We then employed the stacking model with CBR as the meta-model to generate seasonal drought maps based on SPEI-3 across different years and seasons. Our predicted SPEI-3 was compared with drought maps generated by the GPCC and ERA5 reanalysis datasets. We found a strong correlation between them with R2 exceeding 0.8 across multiple years and seasons, indicating our machine learning stacking ensemble approach had a good performance in monitoring seasonal drought conditions. The findings of our study establish a standardized protocol for predicting droughts in the agricultural sector, enabling its widespread application irrespective of variations in terrain and climate.